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The quark-hadron duality in CC and NC neutrino interactions is discussed under assumptions 
that single pion production is described accurately by the Rein-Sehgal model and that it allows 
reconstruction of the inclnsive cross section in the resonance region. The duality is measured by 
means of integrals of structnre functions in the Nachtmann variable for < 3 GeV^. The results 
depend on the precision with which contributions from single pion production channels in the overall 
cross sections are known. Several approaches to evaluate them are compared. The duality is 
predicted to be seen for proton target reactions and to be absent for neutron and isoscalar targets. 
Two-component duality between resonant and valence quark contributions to structure functions is 
also investigated. 
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1. INTRODUCTION 

In recent years the quark-hadron (QH) duality has been a subject of many experimental and theoretical studies. 
The idea of duality comes from Bloom and Gilman Q . It was based on the comparison of eN F 2 resonance structure 
functions at small values of with their DIS counterparts at large Q^. Bloom and Gilman observed that in plots in 
the variable {2Mv-\-M "^)/resonance peaks are approximately averaged by the DIS structure function. More recent 
experimental data on duality were reported in Q where several plots of F 2 for in the range from 0.3 to 3.3 GeV^ 
as functions of the Nachtmann variable ^ were presented on the same figure. The resonance peaks for ^ > 0.2 are 
seen to be averaged by the DIS structure function calculated at = 5 GeV^. The QH duality can be analyzed in 
quantitative way by calculating the ratios of integrated strengths over a range in ^ covering a chosen set of resonances 
and the DIS structure function. The agreement is on the level of 10%. 

In the theoretical analysis of the duality Q one introduces a language of twist expansion of moments of structure 
functions in powers of 1/Q^ . The duality then means the suppression of higher twists Q. A possible explanation 
how the coherent amplitude in the resonance region can be equal to incoherent sum of amplitudes from the quark 
constituents of nucleon is proposed in in terms of SU(6) symmetry by means of cancellation of contributions of 
positive and negative parity. The knowledge of relative strengths of electromagnetic and neutrino induced N —> N* 
transitions leads to theoretical predictions concerning the domain in W in which the duality should hold. The 
consequences of SU(6) breakdown for ratios of unpolarized and polarized structure functions are also investigated in 

The QH duality is usually discussed theoretically and analyzed experimentally in the context of eN interactions. But 
the subject is relevant also for neutrino physics. Here the exact data is missing and arguments based on duality can be 
used to provide better estimates of ^N structure functions and cross sections in the few GeV neutrino energy region 
where they are known with insufficient accuracy Q- For these energies it is necessary to consider both quasi-elastic 
and inelastic channels. Single pion production (SPP) channels are typically treated separately from more inelastic 
ones which are accounted for by extrapolating DIS formalism as much as necessary. This approach carries a lot of 
uncertainty and requires a better theoretical understanding in the whole kinematical region. Such understanding and 
control over the numerical procedures can follow from the QH duality analysis. 

When trying to discuss the QH duality in the neutrino interactions the main obstacle is a lack of precise experimental 
results. The existing SPP data is poor in precision and statistics [^. In the future the data will hopefully become 
precise enough to impose more rigid constraints on theoretical models but for a moment the natural strategy 
is to analyze in detail a generally accepted model. From the point of view of Monte Garlo application this is the 
Rein-Sehgal (RS) ^ model. 
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The RS model was developed to describe SPP induced by le^on-nucleon NC and CC interactions. It is currently 
used in almost all neutrino Monte Carlo generators of events 13- It is based on the quark model computations of 
the hadronic current as proposed and developed by Feynman, Kislinger and Ravndal |ll| . In the original model the 
contributions from 18 resonances with masses smaller than 2 GeV are added in the coherent way. The scattering 
amplitudes contain several form factors. The vector form factors are obtained from the electro-production data by 
means of standard CVC argument. Axial form factors are less constrained by theoretical arguments (PCAC) and the 
available data. The RS model contains also an ad hoc prescription for computing the non-resonant background. The 
extra amplitude is introduced with the quantum numbers of the Pn resonance but without the Breit-Wigner term and 
then added in the incoherent way. Its strength is fine tuned in order to obtain an agreement with the experimental 
SPP data. The RS model can be supplemented by containing terms na (m is the charged lepton mass) absent 
in the original paper. 

The discussion of the QH duality for neutrino scattering based on the RS model is in a very important point 
different from the electron scattering analysis based on the experimental data. The aim of the RS model is to describe 
SPP channels only. In the electron scattering studies the experimental data is that of the inclusive cross section. In 
order to perform the analogous analysis for the RS model it is necessary to extract and add contributions from more 
inelastic channels in the kinematical region of invariant hadronic mass up to 2 GeV. 

Our analysis is based on two basic assumptions. The hrst is that the RS model predictions for SPP cross sections 
are fairly close to what will come out from future precise cross section measurements. The second is that we know 
the probability that at a given point in the kinematically allowed region the final state is that of SPP. Throughout 
this paper we will call the probabilities 1-pion functions. We obtained these functions numerically M using the 
LUND algorithm ^3. The functions (one for each exclusive SPP channel) were discussed also in they turn 

out to depend only on the invariant hadronic mass. For the value of invariant mass kF ~ 1.6 GeV the contribution 
from more inelastic channels amounts to about 50%. Using these functions one can rescale SPP contributions to the 
structure functions and evaluate the overall structure functions in the region lU < 2 GeV. In section IT^ we discuss 
the precision with which the 1-pion functions are reconstructed. We present the experimental data for the proton 
photoproduction 1 -pion function 113 We present also an alternative approach to calculate the 1-pion function for 
available hadron multiplicity data and KNO model |l^ . 

The aim of our paper is to investigate the question if the duality holds also in z/N reactions. Our methodology is to 
repeat the analysis done for eN scattering. We calculate structure functions as they are defined by the RS model and 
compare with the DIS structure functions based on GRV94 PDFs and evaluated at Qjjjs — GeV^. We perform 
also a quantitative comparison of integrated strengths over resonances. All the comparisons are done first for proton 
and neutron structure function and then for their average i.e. for isoscalar (deuterium-like) target. Both charged 
current (GG) and neutral current (NG) reactions are discussed. 

It is well known that GG SPP channels on proton and neutron have distinct properties. The strength od A 
resonance for proton reaction is three times as big as its neutron reaction counterpart due to isospin rules. In the 
neutrino-proton reaction there is no need to introduce non-resonant background which gives significant contribution 
to the neutrino-neutron SPP channels. The way in which non-resonant background is treated in the RS model 
in not completely satisfactory. For this reason we address the idea of two-component duality proposed by Harari 
and Freund 0 . They suggested that resonance and non-resonance contributions to the low energy ttN scattering 
amplitude (s-channel) correspond to contributions given by the high energy amplitudes (t-channel) due to Reggeon 
and Pomeron exchange respectively. Using the modern language it can be expressed as existence of a relation between 
resonance/valence quark and non-resonant/sea quark contributions to the structure functions. A confirmation of this 
idea in eN interactions was found in 0 : the F 2 structure function averaged over resonances at low values of the 
Nachtmann variable (^ < 0.3) behaves in the way which strongly resembles the behavior of valence quark contribution 
to DIS scaling curve. There is also a striking similarity between the above mentioned averaged F 2 rescaled by a factor 
of ^ and the xF^ z/N data. If resonance contribution to structure function is dual to the valence DIS contribution 
and if the overall duality is satished then also non-resonant background should be dual to the sea quark contribution. 
The last duality could be then used to provide a model for non-resonant background. 

The question of the QH duality in z/N interactions was discussed already by Matsui, Sato and Lee 0 | and by 
Lalakulich, Paschos and Piranishvili 0|. In 0| the Lee-Sato model [0 for A production in eN and z/N scattering 
was analyzed. It was shown that in the vicinity of A excitation peak the local duality holds for both proton and 
iso-scalar structure functions in GG and NG neutrino reactions. The resonance model investigated in |2l| includes four 
resonances P 33 , Pn, D 13 and S'!! of lU < 1.6 GeV. The model (unlike the RS model) contains correction for the non¬ 
zero charged lepton mass. The conclusions are in agreement with those contained in [0. Some qualitative elements 
of the present analysis of the RS model can be found in [0. In several theoretical and practical issues related 
to the problem of how to combine smoothly the RS and DIS contributions in Monte Garlo generators are addressed. 
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FIG. 1: Elasticities of resonances 


and the 1-pion functions for neutrino CC reactions. 
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FIG. 2: Comparison of the 1-pion functions for neutrino-proton CC scattering (solid line) and proton photoproduction (dashed 
line). The 1-pion function for the photoproduction was extracted based on the data from 0 


Using the idea of duality one can combine experimental and theoretical arguments and find suitable modification of 
the structure functions which average over resonances for small values of Q^. One expects that duality should hold for 
> 0.5 GeV^ 0. For approaching zero F 2 structure function for electro-production should behave like F 2 ~ 
due to gauge invariance. The presence of axial current modifies this behavior. An important piece of information is 
provided by the Adler sum rule [ 2 ^. It is argued that vector and axial parts of the structure functions should be 
modified in a different way . 

Our paper is organized as follows. In Section 2 the necessary theoretical introduction is given. The recipe to obtain 
structure functions from the Rein-Sehgal model is presented and the 1-pion functions are introduced. The functions 
7 ?.2,3 which measure how well the duality holds are also defined and ambiguities in their definition are discussed. In 
Section 3 we present the results of the numerical analysis and their discussion. Section 4 contains the conclusions. 


2. THEORETICAL FOUNDATIONS OF THE NUMERICAL ANALYSIS 
2.1. Rein-Sehgal structure functions 

We consider the following SPP charged current and neutral current reactions: 

i^{k) +J\f{p) l{k') I + N' + TT 

v{k) + N{p) ^{k') + N* {p) ^ V + N' + TT 

In the LAB frame the momentum transfer is: 


= V 


2 


= k^-k'^ = q^q^ 


q^ = -Q\ 


( 1 ) 
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FIG. 3: Dependence of the Nachtmann variable ^ on hadronic invariant mass calculated at Q = 0.4, 1, 3 and 10 GeV^ 


The leptonic current is defined as: 

^teptou = - -i5)u(k). 

In the RS model the leptonic mass is set to be zero. In this limit 

^P-^lepton 

One can introduce the basis of three vectors of length ±1 orthogonal to q^: 

1 


el = —(0,l,-i,0), 
.A. _ 2 

V 2 

1 


= —(0>-l>-b0), 


60 — 


Correspondingly, the leptonic tensor can be decomposed as: 


(g,0,0,^). 


- g^^'k ■ k' - = 


( 2 ) 

(3) 


(4) 


= ^ (5) 

a,l3e{S,L,R) 

When we calculate the contraction of the leptonic tensor with the hadronic tensor 



Ws), 

(6) 

{M is the nucleon mass) we find that 

L^’'W,, = LZ,W,,, 

(7) 

where 



rpi' 

^ diag 

= A^el{eir + B^el{eir + C^el{e‘ir. 

(8) 

^ B^, are Lorentz scalars which 

can be evaluated in the LAB frame: 


A" 

= L^.elieir = ^ {{2E - rY - q^) , 

(9) 


= L^^elieir = ^{2E-R + qY, 

(10) 


= L,^el{eir = ^{2E-R-qY. 

(11) 
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Decomposition of the leptonic tensor entails the decomposition of the cross section into three contributions ctl, an, as 
which are interpreted as cross sections of intermediate boson in given polarization states scattered off nucleon. In its 
final form the Rein-Sehgal formula for the cross section reads: 


where 


Thus we can write: 


Gp cos^ Oc 
47r^ 



+ 2uvas) dQ^dv, 


2 Ij2 2 pi2 9^ /12 


^ ^ + ^^crs) dQ'^dv. 


( 12 ) 


(13) 


'^L,R,s are then calculated within the quark model and are given in the explicit way for each SPP channel separately. 
In order to identify aL,R,s as linear combinations of the structure functions we calculate: 


(W2^ - IPi) + (iTi + ITg^) + (iTi - 


and 


d\ 


dvdQ'^ 


cos^ 9c 


(ir^ - «'0 + (»'. + "-iw) + ("-■. - 


A simple comparison gives the following identification of the tr's in terms of IT's: 

„2 


- ITi + IT 2 




= -CTS, 




(14) 


(15) 

(16) 
(17) 


Using the definition of the structure functions F's we write down the final expressions for the RS model structure 
functions: 


fRS 

= MWi = M^{aLEaR), 

ZTT 

(18) 

pRS 

K, Q‘^ 

= iyW 2 = v- -^ (2cts + ctl + an), 

ZTT 

(19) 

pRS 

= vWz = v- — {aL - an). 

TT q 

(20) 


2.2. DIS structure functions 

We apply the simple model for the DIS structure functions: they are given by appropriate linear combinations of 
the parton distribution functions (PDFs). In the kinematical region we are interested in the charm contribution can 
be neglected. For CC i/N interaction the DIS structure functions are |2^: 



= 2x (d cos^ 9c + s sin^ 9c + u) , 

(21) 

xF^^{rp) 

= 2x [d cos^ 9c + s sin^ 9c — u) , 

(22) 

F^^{vn) 

= 2 x (u cos^ 9c + s sin^ 9c + d) , 

(23) 

xFf^ii^n) 

= 2x [u cos^ 9c + s sin^ 9c — d) ■ 

(24) 
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For i/N NC interaction the DIS structure functions are: 


— 2a: + glj) (u + w) + (d + d + 2s)^ , (25) 

xF^^{vp) = 2x (^{gl - g\) {u - u) + {g'\ - g'\) {d - d)^ , (26) 

= 2x {{g\-\-g\) (^d + d)-\-(g'+ g'{u + u-\-2s'^ , (27) 

xF^^{vn) = 2x {{gl -g\){d-d) + {g'\ - g'\) {u - m)) . (28) 

with 

12 2 

9l =-^ --^sin^ 6w, gn =9w, (29) 

fl'i = -1 + 1 sin^ Ow, g'R = 1 sin^ Ow (30) 

Wherever we present the plots of Fi we use the Callan-Gross relation: 

F2=2xFi. (31) 

In the quantitative analysis we restrict ourselves to F 2 and xF^ only. 


We use GRV94 (LO) PDF’s which are defined for > 0.23 GeV^ and a; > 10 ® and distinguish valence and 
sea quark contributions. GRV94 PDF’s are used in many Monte Garlo generators of events Q- 


2.3. 1-pion functions 


The 1-pion functions are defined for each SPP channel separately as probabilities that at a given value of W the 
final hadronic state is that of SPP: 


MW) ^ (32) 

dW 


The 1-pion functions used in this paper were obtained from the Monte Garlo simulation based on the LUND algo¬ 
rithm. Therefore they are defined by fragmentation and hadronization routines implemented there. The comparison of 
the 1-pion functions with elasticity factors of resonances included in the RS model is shown in Fig. ^ The agreement 
is satisfactory. It was also checked that the simulations based on LUND give rise to charged hadron multiplicities 
consistent with the experimental data |^. 

We assume that the following relation holds between the structure functions of the RS model and the overall 
structure functions for the inclusive cross section: 


FPPP{x,Q^) 


FfHx,Q^) 

fiAWix,Q^))’ 


(33) 


where j = 1, 2, 3. In the case of neutron (GG reactions) and NG structure functions in the above formula we apply 
the sum of the 1-pion functions for two SPP channels. 

We will see that our results are sensitive to the details of the 1 -pion functions in the region of W G (1.5, 2) GeV 
where the rescaling effects are most important. We investigated this point in more detail: 

i) We extracted the 1-pion function from jp photo-production data [13, see Fig. |21 We conclude that it is rather 
similar to the function we used in our numerical computations. 

ii) We tried to evaluate the 1 -pion functions from available hadron multiplicities data in neutrino reactions |l 6 | ex¬ 
trapolating the predictions to the region W G (2, 3) GeV. 

We know the average charged hadron multiplicities 
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{nch)^p = -0.05 ± 0.11 + (1.43 ± 0.04) \n{W'^), 


{nch)un = -0.2 ± 0.07 + (1.42 ± 0.03) ln{W^), 
and the neutral pion multiplicity measured in vp reactions: 


(n^o)^p = 0.14 ± 0.26 + (0.5 ± 0.08) \n{W^). 

We assume KNO distribution of multiplicities and that in the first approximation there are only nucleons and pions 
in the final state. For vp reaction it is also necessary to make some assumptions about charged pions. We expect 
that the fraction of charged pions is the same as in the pp reaction for which the experimental data is available. 
We obtained the following values for the 1-pion functions at IF = 2 GeV: 0.14 for lyp and 0.38 for i/n. The obtained 
values must be reduced by ~ 10 % due to the presence of other exclusive channels. We conclude that the KNO 
results seems to be in general agreement with basic properties of our 1-pion functions: the fraction of SPP channels 
on neutron is bigger then on the proton SPP and the orders of magnitude are quite similar. 


2.4. Kinematics 


In the QH duality analysis different kinematical regions are simultaneously involved in the discussion. 

A common presentation of the duality is done by means of the comparison of plots of structure functions in the 
Nachtmann variable: 




2x 

1 + y/l + Ax^M-^/Q^ 


(34) 


which takes into account target mass corrections. 

The resonance region is defined in terms of invariant hadronic mass as IF G (M + ttItt, 2 GeV) which is the natural 
choice for the Rein-Sehgal model. Other options {Wmax < 2 GeV) for the definition of the resonance region were 
considered in |^. 

Fig.El illustrates the dependence of the Nachtmann variable on hadronic invariant mass at fixed values of Q^. For 
typical for the resonance region G (0.5, 3) GeV^ one obtains the following domain in the Nachtmann variable: 
^ € (0.13,0.76). This region in ^ when combined with ^0 GeV^ corresponds to: 


O' 

= 5 GeV^ 

IF e (1.3, 5.8) GeV, 

(35) 

o' 

= lOGeV^ 

IF G (1.8, 8.2) GeV, 

(36) 

o' 

= 20GeV^ 

IF G (2.5, 11.6) GeV. 

(37) 


In our numerical analysis we use Q]jjg = 10 GeV^. 


2.5. Quark-Hadron Duality 


The QH duality is said to be present on the quantitative level if the following relation between resonance and scaling 
structure functions holds: 

Ql^s) « QIjs). (38) 

The above equation should hold for different values of Q\es characteristic for the resonance production and for a 
fixed value of Q^jg- The region of integration ~ RES region - is defined to be identical with the resonance region of 
the RS model: Wmin = M + m-,^ and Wmax = 2 GeV, which is then translated into appropriate region in Q 


^min 


^{Wmax, Q'rEs)^ 


^max 


^(Wmin,QREs) 


(39) 





FIG. 4: Uncertainties in due to different definitions of Q\>is- Solid line corresponds to on) and dashed line to 
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FIG. 5: Comparison of the Rein-Sehgal structure functions at = 0.4, 1 and 2 GeV^ with the appropriate scaling functions 
at (5i)/s=10 GeV^. In the first row xFi, F 2 and xFs structure functions for CC neutrino-proton scattering are plotted. In the 
second row the structure functions for CC neutrino-neutron scattering are shown. 


In the quantitative analysis we define ratios of two integrals over the resonance region: 


n{f,Ql;g,Ql) 



dUi^Qn) 


(40) 


^min 


d^9{^,Ql) 



























9 


xFJ 


xFo' 





xf; 


xFo' 





FIG. 6: The same as in Fig.|21but for NC reactions. 


Using the above quantity, we define the function: 

'^2{Qres^Qdis) — ^QreS’tQdis) ■ 


(41) 


and 


'^^{Q'res^Q'dis) — ^Q'res'^^^F^^^Q'dis) ■ 

(42) 

There is an ambiguity in the quantitative analysis of duality because the above functions depend on the arbitrarily 
chosen value of The differences between scaling curves calculated at different are not relevant for making 

qualitative statements about the duality but do matter in quantitative analysis. To illustrate this we calculate 

7 ^l0/5 = n = 10; uf Qhs = 5), (43) 

and 

7 ^l0/20 = F Qhs = 10; F^^^, QIjs = 20), (44) 

the integration region of the above integrals is defined by Q\es- results are shown in Fig. ^ We see that the 
ambiguity is of the order of 7% and is largest for Qrer ~ 1.5 GeV^. This limits the precision of the quantitative 
statements about the QH duality. 

In the investigation of two-component duality we single out resonant and non-resonant contributions to the RS 
model structure functions 




RES 


— Fj^res F Fj 


( 45 ) 
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FIG. 7: Comparison of the Rein-Sehgal structure functions at = 0.4, 1 and 2 GeV^ with the appropriate scaling functions 
at Qi)/s=10 GeV^. In the first row the plots of the xFi, F 2 and xFs structure functions for CC neutrino-isoscalar target 
scattering are presented. In the second row structure functions for NC neutrino-isoscalar target scattering are shown. 


and we also separate valence and sea quark contributions to the DIS structure functions: 

pDis _ p 4_ p 

^j — ^j,sea I 

We calculate the following functions: 

T^2‘^’'iQ'RESyQ'DIs) ='^{P2,res,Q'}iES'^^'^,valyQ‘DIs) ' 

and 

'^3°‘’'{Q'RESyQi)Is) = '^ {^P3,res,Q'}iES'^^^3,val,Q'DIs) ' 


(46) 

(47) 

(48) 


3. NUMERICAL RESULTS AND DISCUSSION 

In the numerical analysis we confine ourselves to the case of neutrino interactions and leave out the antineutrino 
ones. 

In Figs. 0 0 we present a comparison of the scaling structure function with the RS structure functions calculated 
at Q\es ~ 0-4, 1 and 2 GeV^. The Figs. 0 and El correspond to CC and NC reactions respectively with proton 
structure functions in the upper row and neutron structure functions below. 

In the case of the RS model for neutrino-proton CC reaction the A resonance contribution dominates overwhelmingly 
over other resonances. One can see the typical manifestation of local duality: the sliding of the A peaks (calculated 
at different Q\es) along the scaling function. 

For neutrino-neutron CC reaction the resonance structure is much richer. The contributions from the A are usually 
dominant but those from more massive resonances are also significant. In the figure with the structure function 
three peaks of comparable size are seen. The DIS contributions dominate over the RS ones in this case. 
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FIG. 8: The comparison of the proton and neutron CC F 2 structure functions for the RS model (dashed line), rescaled RS 
model (dotted line) and the DIS (solid line). The RS structure functions are calculated at = 1 GeV^ while the DIS curves 
are plotted for = 10 GeV^. 


Our plots are comparable with those obtained in Ref. I^^]. However in our paper the strengths of A peaks for all 
the structure functions are proportionally lower. The difference between structure functions based on CTEQ6 |2n| 
and GRV94 PDF’s is very small. 

It is seen that it is virtually impossible to have simultaneously duality in CC reactions on proton and neutron 
targets. The strength of A excitation on proton is approximately three times as big as for neutron and DIS cross 
section on neutron is much bigger then on the proton. 

For the NC reactions the dominant A peaks also slide along scaling function. The structure functions for the proton 
and neutron are almost the same. 

In Fig.[7|the analysis is performed for the isoscalar target. The plots of the structure functions for CC (upper row) 
and NC (lower row) reactions are presented. In each plot the A peak slides along the scaling function and the local 
duality is seen. 

Apparently (with the exception for CC reaction on proton) there is little hope for the QH duality in the whole 
resonance region: the scaling structure functions are on average larger then the RS ones. Only the local duality is 
present after a suitable region in W around Ma = 1.23 CeV is chosen. But in the Figs. 0- El the rescaling of 
RS structure functions by means of the 1-pion functions (see Fq. CT ) is not yet included. The rescaling procedure 
increases the RFS structure functions making the duality more likely to appear. We notice also that the rescaling 
cannot spoil the statements about the local duality around the A resonance because the values of the 1-pion function 
for W in the vicinity of Ma are close to 1 (see Fig. P). 

In Fig. 1^ we show how the resonance structure functions are modified by means of the 1-pion functions. The 
modifications apply mainly to hadronic invariant mass close to 2 CeV. 

In order to perform a quantitative analysis of the duality we make use of the functions TZi defined in Fqs. 141I48II . 
We restrict our plots to the values of ^ GeV^ characteristic for the resonance production. 

In what follows we use the RFS structure functions rescaled by means of the 1-pion functions. We have checked that 
introduction of the 1-pion function improves the duality significantly. For example for CC reaction on neutron TZ 2 is 
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Charged Current 



0 0.5 1 1.5 2 2.5 3 

[GeV^] 


FIG. 9: The functions TZ 2 for different targets and reactions. The ratios are calculated for CC and NC structure functions in 
the cases of proton (solid lines), neutron (dotted lines) and isoscalar target (dashed lines). 


increased by a factor of ~ 1.55 and for proton by ~ 1.39. The difference is caused by the overwhelming dominance 
of the A excitation in the case of proton. 

A characteristic feature of most of the plots of T^jiQjiEs) ^ presence of two qualitatively distinct behaviors. For 
Qj^ES smaller then ~ 0.5 GeV^ the functions TZj vary quickly while for larger values of Q\es they become slowly 
changing. This seems to correspond to predictions done in ra. Our statements about the duality will apply only to 
the region of > 0.5 GeV^. 

In Figs. El and El the plots of 1^2 and for proton, neutron and isoscalar targets are presented. In the case of GG 
interaction the duality is seen on the proton target (accuracy < 20%) but for the neutron and isoscalar targets the 
duality is absent. In both cases the average strength of resonance structure functions amounts to only about a half 
of the strength of DIS structure functions. The plots for the NC interactions are almost independent on the target 
and in all the cases the DIS contributions are approximately two times as big as resonance ones. A different choice of 
Qdis^ namely Qjjjg = 20 GeV^ makes the values of 7?.2,3 even lower (see Fig. 

The remaining plots address the question of two component duality. We concentrate on the case of the possible 
duality between the resonance and valence quark contributions. 

In Fig. m the plot of 7?.™^ for the CC interactions is shown. We notice the good duality picture in the case of 
proton target but a huge departure from duality in the case of neutron and isoscalar targets. It is worth noting that 
this discrepancy is larger than one shown in Fig. El where the general (not two component) notion of duality was 
discussed. The novel feature is the apparently singular behavior at low Q\es'- rises quickly in contrast with TZ 2 
falling down when Q\es approaches zero. 

The explanation of this follows from the Fig. II21 where the region of small Q^es analyzed in more detail. We 
notice that for Q'j^^ES approaching zero the valence quarks scaling function tends to zero while the resonance strengths 
remains virtually unchanged. 

Finally in Fig. El the analogous two-component duality analysis is done for The discussion of xF^ seems to 

be favorable for the two-component duality because in the DIS contribution on the isoscalar target there is no sea 
quark contribution. We remind also that for the GG reaction on the proton the non-resonant contribution is absent. 
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Charged Current 



0 0.5 1 1.5 2 2.5 3 
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FIG. 10: The same as in Fig. |^but for xF^ (ratio TZs). 


Charged Current, resonant to valence 



FIG. 11: The plots of functions 7?.2“* defined in Eq. The computations are performed for the GC reactions for proton (solid 
line) neutron (dotted line) and isoscalar targets (dashed line). 


In Fig. El we see that two component duality is satisfied within ^30% for the proton target but it is absent for 
neutron and isoscalar targets. We notice also that contrary to what we have seen in the plots for T?.™* now at low 
Q\es curves tend to zero. 

The explanation of this behavior follows from the Fig. El One can see that in the case of xF^ both the resonance 
and valence quark structure functions fall down for approaching zero. The behavior of xF^ is the same as that 
discussed in [l^ . 

We do not present plots exploring the duality between the non-resonant part of the resonance model and the sea 
quark contribution. No sign of two component duality is seen in this case. 
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F P 
^2 




FIG. 12: Comparison of the proton and neutron CC Rein-Sehgal structure functions at = 0.05, 0.1, 0.2 and 0.4 GeV^ 
with DIS structure (dashed line) at Q|j/s=10 GeV^ and the contribution valence quarks (solid line). 


Charged Current, resonant to valence 



FIG. 13: The same as in Fig. I1 1 1 hut for xFs. 


4. CONCLUSIONS 

The quark-hadron duality in the vN reactions has been investigated by comparing the structure functions obtained 
from the Rein-Sehgal model and those from the deep inelastic formalism. The 1-pion functions were used to construct 
comparable quantities. The qualitative analysis was based on the plots of the RS and the DIS structure functions 
at several values of while the quantitative one was based on the functions TZj defined as the ratios of the 

corresponding integrals. 

We are aware that our model of resonance structure functions is a subject of several uncertainties. Here is a list of 
them: 
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Isoscalar target: xFj 



FIG. 14: Comparison of the isoscalar CC Rein-Sehgal xFs structure functions at — 0.05, 0.1, 0.2 and 0.4 GeV^ with 
appropriate valence quark contribution to scaling function (solid line) calculated at Q\,js = 10 GeV^. 


ratios for F2 (CC), isoscalar target 



ratios for F2 (NC), isoscalar target 



FIG. 15: The ratios presented in Fig. I^for isoscalar target are calculated without 1-pion functions but with all the elasticities 
equal to 1 with amplitudes added in either coherent or incoherent way. 


a) It is possible that the structure functions extracted from the RS model are underestimated for VF > 1.7 GeV where 
tails of heavier resonances are not included . 

b) In the original RS model the non-resonant background is treated in not a satisfactory way. In particular the shape 
of the non-resonant background does not seem to agree with the precise electro-production fits presented in Ref. 

ga. 

c) Our conclusions depend on the precision with which we reconstructed 1-pion functions. Therefore we decided also 
to investigate the duality with no assumptions about the 1-pion functions but rather with overall cross sections for 
the resonance production (i.e. under the assumption that all the elasticities are equal 1). We considered two cases: 
(i) the resonances are added incoherently and (ii) the interference patterns are the same as in SPP channels. The 
typical results are shown in Fig. 1151 predictions of all three models turn out to be similar. 
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Our main conclusions can be summarized as follows: 

1) for > 0.5 GeV^ the duality is present in the whole resonance region W G {M + to^, 2 GeV) with an accuracy 
of ^ 20% only for GG proton target reaction. We remind that the way in which the duality is defined carries an 
uncertainty of ~ 5%; 

2) from the Figs. El-[3 it follows that there is also a local duality for isoscalar target in the case of CC reaction and 
for all the targets in the case of NG reactions in a suitable chosen vicinity of the A resonance. 

The results obtained in this paper can be useful for the investigation of the question how to modify DIS structure 
functions in the low resonance region so that they provide a good average description |82|| . Such modifications 
should be confronted with available neutrino scattering data from CHORUS, NOMAD and NuTeV experiments [^ . 
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